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Abstract 

The nucleon strange vector and isoscalar electromagnetic form factors are studied us- 
ing a spectral decomposition. The KK contribution to the electric and magnetic radii as 
well as the magnetic moment is evaluated to all orders in the strong interaction using an 
analytic continuation of experimental KN scattering amplitudes and bounds from unitar- 
ity. The relationship between non-resonant and resonant KK contributions to the form 
factors is demonstrated, and values for the vector and tensor <f>NN couplings are derived. 
The KK spectral functions are used to evaluate the credibility of model calculations for 
the strange quark vector current form factors. 
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1 Introduction 



The reasons for the success of the constituent quark model of light hadrons remains one of 
the on-going mysteries of strong interaction physics. Although deep inelastic scattering has 
provided incontrovertible evidence for the existence of gluons and QCD current quarks in the 
lightest hadrons, these degrees or freedom are manifestly absent from the the quark model. 
Nevertheless, a description of light hadrons solely in terms of constituent quarks moving in 
an effective potential has been enormously successful in accounting for the mass spectrum and 
other properties of low-lying hadrons. Various explanations for this situation have appeared 
in the literature, including the simple and intuitive idea that the sea quarks and gluons of 
QCD "renormalize" the valence current quarks into the constituent quarks of the quark model 
IJ. In this picture, for example, the multitude of QCD degrees of freedom appear to a long 
wavelength probe primarily as single objects carrying the quantum numbers and effective mass 
of the constituent quark. From the standpoint of the quark-quark effective potential, gluon 
and sea-quarks are similarly un-discernible - as they help renormalize the quark model string 
tension into the physical value used as model input ||. In fact, most low-energy observables 
studied to date are unable to uncover explicit signatures of QCD degrees of freedom. 

There have been, however, a few exceptions to this situation. Of particular interest are 
observables sensitive to the presence of strange quarks in the nucleon. In contrast to up- 
and down-quarks, which appear both as valence and sea quarks, strange quarks constitute a 
purely sea-quark degree of freedom. Being the lightest such objects, they ought to generate the 
largest effects (in comparison to the heavier quarks). Consequently, nucleon matrix elements of 
strange quark operators provide an interesting window on the qq sea and as such may shed new 
light on the connection between non-perturbative QCD and the quark model. Indeed, were 
strange-quark observables found to be vanishingly small, one might ascribe the quark model's 
success partly to the numerical insignificance of sea quark effects. [] In fact, the situation is 
more ambiguous. As is well-known, analyses of the "cr-term" in nN scattering, the g\ sum in 
polarized deep inelastic scattering, and N deep inelastic scattering suggest that non- 

trivial fractions of the nucleon mass, spin, and light-cone momentum arise from the ss sea (see 
Ref. |P] and references therein). Evidently, the most naive explanation for the quark model's 
validity is ruled out by these analyses. 

More recently, a well-defined program has begun to determine the matrix element (Nlsj^slN) 
using parity-violating elastic and quasi-elastic electron scattering from the proton and nuclei 
H]. The first result for the magnetic form factor associated with this matrix element has been 
reported by the SAMPLE collaboration at MIT-Bates [|]: 

G s M (q 2 = -0.1 GeV 2 ) = 0.23 ± 0.37 ± 0.15 ± 0.19 , (1) 

where q 2 is the four-momentum transfer squared. (The first error is statistical, the second is 
the estimated systematic error, and the last uncertainty is due to radiative corrections entering 

3 The reason why non-perturbative QCD produces small sea-quark effects at low-energies would remain to be 
explained, however. 
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the analysis.) Although the value is consistent with zero, the error bars are large. Improved 
accuracy is expected when the full data set is analyzed. Similarly, a combination of the strange 
magnetic and electric form factors have been determined by the HAPPEX collaboration 0: 

G S E + 0.39G s M (q 2 = -0.48 GeV 2 ) = 0.023 ± 0.034 ± 0.022 ± 0.026 , (2) 

where the first two errors are again of statistical and systematic origin, respectively, and the 
last one arises from the estimated uncertainty in the electric neutron form factor. While no 
definitive conclusion can as yet be made regarding the experimental scale of {N\s^ffj,s\N), one 
expects to be able to do so at the conclusion of the measurements. 

In contrast, the theoretical understanding of (iV|s7 M s|iV) is much less clear. The difficulty 
lies with the mass scales relevant to strange quark dynamics. In contrast to the heavy quarks, 
for which m q ^> Aq C d, the strange quark has m s ~ Aq C d- Consequently, the lifetime of a 
virtual ss pair is commensurate with typical strong interaction timescales, allowing the pair to 
exchange a plethora of gluons with other quarks and gluons in its environment. The dynamics 
of the pair are therefore inherently non-perturbative. Given the present state of QCD theory, 
a complete, first principles treatment of (Nlsj^slN) has remained beyond reach. Attempts to 
obtain this matrix element on the lattice have produced two results for G s M (q 2 = 0), neither of 
which agree with each other nor with the first SAMPLE results [B, |7|, and one result for the 
slope of G% at the origin with large error bars [0]. 

An alternative - and more popular approach - has been to employ various effective frame- 
works, with varying degrees of model-dependence. These frameworks have included nucleon 
models, chiral perturbation theory (ChPT), and dispersion relations. Generally speaking, the 
degrees of freedom adopted in each of these approaches have been hadronic rather than quark 
and gluon, given that the lifetime of an ss pair permits it to form strange hadronic states. 
Apart from a few exceptions, effective approaches do not address the way in which QCD sea 
quarks hadronize. Hence, the connection with QCD is indirect at best, with each approach 
emphasizing some aspects of the strong interaction to the exclusion of others. 

Not surprisingly, the range of predictions for the strangeness form factors is broad. In 
particular, the breadth of model predictions appears to be as wide as the variety of models 
that has been used even though the same models are in reasonable agreement for standard 
nucleon observables 0, ^ [|. This situation illustrates the sensitivity of sea quark observables 
to model assumptions and the limited usefulness of models in making airtight predictions. One 
might have hoped for more insight from ChPT, which relies on the chiral symmetry of QCD 
to successfully account for a wide variety of other low-energy observables |IJ. Unfortunately, 



ChPT is unable to make a prediction for the leading non-vanishing parts of G S M or G S E since 



the leading moments depend on unknown counterterms JTTJ . Recently, however, it has been 
noticed that slope of G S M at the origin is independent of unknown counterterms to 0(p 3 ) |i~2|j . 

In the present study, we turn to dispersion relations (DR's) to derive insight into (Nls'j^slN). 
Like ChPT, DR's rely on certain general features of QCD (and other field theories) to relate 
existing experimental data to the observables of interest. In the case of DR's it is analyticity and 
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causality, rather than chiral symmetry, which allow one to make the connection. Although DR's 
do not bear on the way in which QCD quarks and gluons form intermediate strange hadronic 
states, they do provide an essentially model-independent framework for treating the way in 
which those states contribute to the form factors. We view them as providing an intermediate 
step toward understanding the strange quark form factors at the fundamental level of QCD. 
Because of their generality, they also allow us to evaluate the credibility of several model 
predictions. 

Our use of DR's to study nucleon form factors is not new. The spectral content of the 
nucleon isovector form factors has been clearly delineated using this approach |13, 14]. It is 



now known that both an un-correlated tttt continuum as well as the tttt —>■ p resonance play 
important roles in the low-g 2 behavior of these form factors. The g 2 -dependence of the isoscalar 
EM form factors has been successfully reproduced using DR's under the assumption of vector 
meson dominance (VMD). The results have been used to infer relations between the uNN and 
<j)NN coupling strengths and to make predictions for the strange quark form factors [ITR ITB, |T7 



However, the relationship between the resonance and continuum contributions to these form 
factors has not been previously established. Consequently, a number of model predictions have 
appeared which rely on the assumption that the uncorrelated continuum ( "meson cloud" ) gives 
the largest effect. These meson cloud calculations have generally entailed a truncation at second 
order in the strong hadronic coupling, g - a practice of questionable validity. The corresponding 
predictions have generally been in disagreement with those obtained using VMD. 

In what follows, we consider both the strange quark and isoscalar EM form factors without 
relying on the a priori assumption of vector meson or meson cloud dominance. We focus 
in particular on the contribution from the KK intermediate state. The rationale for this 
focus is twofold. First, the KK state constitutes the lightest intermediate state containing 
valence s and s quarks. Its contribution to the strange quark form factors has correspondingly 
been emphasized in both models and ChPT. Second, given the present availability of strong 
interaction and EM data, the KK contribution can be computed to all orders in g using a 
minimum of assumptions. From an analysis of KN — > KN and e + e~ — > KK data, we show 
that the scale of the KK contribution depends critically on effects going beyond 0(g 2 ) and 
argue that a similar situation holds for the remainder of the form factor spectral content. We 
also 

(a) illustrate the relation between the continuum and resonance contributions, 

(b) evaluate the credibility of several model predictions as well as the 0(p 3 ) prediction of ChPT 
for the magnetic radius, 

(c) derive values for the vector and tensor <f)NN couplings and compare with those obtained 
from isoscalar EM form factors under the assumption of VMD. 

In Refs. ]T8|, [T9|| , we reported on the results of our DR analysis of the KK contribution to the 



nucleon "strangeness radius" (the slope of G S E at the photon point). Here, we expand on that 
analysis to consider the full g 2 -dependence in both the isoscalar EM and strangeness channels 
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and to discuss both the electric and magnetic form factors. Since the DR approach requires 
knowledge of the KK — > NN amplitudes in the unphysical region, some form of analytic 
continuation is needed to complete the analysis. Using backward dispersion relations, we obtain 
the unphysical amplitudes from KN phase shift analyses. The results of this continuation and 
their implications for nucleon form factors constitute a central theme of this paper. 

Our discussion of these issues is organized as follows. After outlining our formalism, we 
perform the spectral decomposition of the form factors and write down DR's in Section |2]. In 
Section [3], we express the spectral functions in terms of KK — > NN partial waves and give 
the corresponding unitarity bounds valid in the physical region of the dispersion integrals. The 
analytic continuation of KN amplitudes which is used in the unphysical region is performed in 
Section f|. A brief description of the analytic continuation and our treatment of the inherent 
problems is given. The results are applied to the nucleon's strange and isoscalar EM form 
factors in Section |. In Section ||, we discuss the contribution of other intermediate states and 
summarize our conclusions. 

2 Spectral Decomposition and Dispersion Relations 

The vector current form factors of the nucleon, Fi(t) and -F^), are defined by 
(N(p')\j,(0)\N(p))=u(p') 

where t = q 2 = (p' — p) 2 . We consider two cases for j^: (i) the strange vector current s^^s and 
(ii) the isoscalar EM current j^ =a) ■ Since the nucleon carries no net strangeness, F( must vanish 
at zero momentum transfer, (i.e. F{ (0) = 0), whereas F 1 (/=0) is normalized to the isoscalar EM 
charge of the nucleon, F± (0) = 1/2. We also define the electric and magnetic Sachs form 
factors, which may be interpreted as the fourier transforms of the charge and magnetic moment 
distributions in the Breit frame, 

G E = F 1 - tF 2 , G M = F, + F 2 , (4) 

with r = —t/4m 2 N . In the case of the strange form factors we are particularly interested in their 
leading moments, the strange magnetic moment and the strange radii: 



Fi(th, 



(p -vY 



u(p) 



(3) 



K 



(r 2 )l = 6 



F°(0) = GtM=V s , 
dGt(t) 



(5) 



dt 



where i = E, M, respectively. A dimensionless version of the radii can be defined by 

2 



Pe 



dG%{r) 



dr 



r=0 



\m 2 N {r 2 y E 



(6) 
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and similarly for p s M . The leading moments of the EM form factors are defined analogously. 

It is conventional to employ a once subtracted DR for F±. Typically, one wishes to predict 
F 2 (0) as well as its t 7^ behavior. In this case, an unsubtracted DR is appropriate. We follow 
this ansatz in the present study and obtain 



F*(t) 



1 



7T Jt 



and 



t r c 
™ = 1 

7T Ju 



ImF£(t') 
t'-t 



ImFf(t') 



dt' 



d£ . 



(7) 



t'(t' - 1) 

As a consequence, subtracted DR's can be written for the Sachs form factors as well. 

From Eqs. it is clear that the quantities of interest are the imaginary parts of the form 

factors. The success of the DR analysis relies on a decomposition of the ImF, into scattering 
amplitudes involving physical states. To obtain this spectral decomposition, we follow the 
treatment of Refs. 



20, 21, 18 and consider the crossed matrix element 



J, = (N(p)N(p)\j,(0)\0) = u(p) 



iF 2 {t) 
2m, 



-a 



{p+pT 



v(p) 



(9) 



where t is now timelike. Using the LSZ reduction formalism and inserting a complete set of 
intermediate states, Im J M may be expressed as 



Im J,, 



71 



:2nf/wY,(N(p)\j N m*)(Mj»mo)v(p)5\p+p-p x ) 



(10) 



where M is a spinor normalization factor and Jat(0) a nucleon source. Eq. ( |10|) determines the 
singularity structure of the form factors and relates their imaginary parts to on-shell matrix 
elements for other processes. The form factors have multiple cuts on the positive real t-axis. 
The invariant mass-squared M\ of the lightest state appearing in the sum defines the beginning 
of the first cut and the lower limit in the dispersion integrals: M\ = t\. Since Eq. (|iy) is linear, 
the contributions of different |A) can be treated separately. 

There is an infinite number of contributing intermediate states |A) which are restricted 
by the quantum numbers of the currents sj^s and j^ =0) [I G (J PC ) = 0~(1 )]. Naively, the 
lightest states generate the most important contributions to the leading moments of the current. 
Moreover, because of the source Jjv(0), the intermediate states must have zero baryon number. 
The lowest allowed states together with their thresholds are collected in Table [TJ Resonances, 
such as the u, do not correspond to asymptotic states and are already included in the continuum 
contributions, such as that from the 3tt state. 

When sufficient data exist, experimental information may be used to determine the matrix 
elements appearing in Eq. (0). However, when the threshold t\ of the intermediate state 
I A) is below the two-nucleon threshold, the values of the matrix element (N(p)\ Jjv(0)|A) v (p) 
are also required in the unphysical region t\ < t < 4m^. In this case, the amplitude must 
be analytically continued from the physical to the unphysical regime. The first cut in the 
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mesonic states 


tA[GeV^] 


baryonic states 


tA[GeV^] 


O/l 


n 1 r 

u. 10 


MM 


o.oo 


5n 


0.49 


N Nnn 


4.67 


7tt 


0.96 


AA 


4.84 


KK 


0.98 


SS 


5.76 


KKtx 


1.28 


AEtt 


5.95 



Table 1: Lowest mass intermediate states contributing to Eq. ([TOD • 



complex t-plane appears at the 3tt production threshold, t = 9m 2 , and higher- mass intermediate 
states generate additional cuts. For example for |A) = \KK) the cut runs from t = Am 2 K to 
infinity. Therefore, the matrix element for KK — > NN is also needed in the unphysical region 
4m 2 K < t < 4:in 2 N , which requires an analytic continuation. 

Some of the predictions for the F? reported in the literature are based on approximations to 
the spectral functions appearing in Eqs. (|7|, §). The work of Refs. [15, [IB], [lTj employed a VMD 
approximation, which amounts to writing the spectral function as ImF?(t) = nYlj cijS(t — m 2 ), 
where "j" denotes a particular vector meson resonance (e.g. uo or 0) and the sum runs 
over a finite number of resonances. In terms of the formalism from above this approxima- 
tion omits any explicit multi- meson intermediate states \n) and assumes that the products 
(iV(p)| J7v(0)|A)(A|s7 M s|0)v(p) are strongly peaked near the vector meson masses. The same 



has been made conventional analyses of the isoscalar EM form factors [13, 14]. 

In contrast, a variety of hadronic effective theory and model calculations for the strange form 
factors have focused on contributions from the two-kaon intermediate state IBL |8|, 01. Even 



though \KK) is not the lightest state appearing in Table [1], it is the lightest state containing 
valence strange quarks. The rationale for focusing on the KK contribution is based primarily on 
the intuition that such states ought to give larger contributions to the matrix element (A^^lO) 
than purely pionic states with no valence s or s quarks. In other words, the kaons represent 
the lightest contribution favored by the OZI rule. Typically, kaon-cloud predictions have been 
computed to 0(g 2 ) only. The results for for p s E in particular are smaller in magnitude than the 
vector meson dominance predictions and have the opposite sign. In what follows, we illustrate 
how both the structure and magnitude of the full kaon cloud contribution differ from the 0(g 2 ) 
result and how a ^-resonance structure appears in the all-orders analysis. 

Although we consider here only the KK intermediate state, we note in passing that the 
validity of this so called "kaon cloud dominance" ansatz is open to question for a variety 
reasons. As can be seen from Table [l], for example, the three-pion threshold is significantly 
below the KK threshold. Consequently, the 3n contribution is weighted more strongly in 
the dispersion integral than the KK contribution because of the denominators in Eqs. (|7|, 
H). Moreover, three pions can resonate into a state having the same quantum numbers as 



6 



the (f> (nearly pure ss), and thereby generate a non-trivial contribution to the current matrix 
element. Indeed, the has roughly a 15% branch to multi-pion final states (largely via a pir 
resonance). Although such resonances do not appear explicitly in the sum over the states in Eq. 
(|lOl), their impact nevertheless enters via the current matrix element (37r|s7 M s|0) and the NN 
production amplitude (N(p)\ Jn(0)\3it)v(p). Thus, the 37r state could contribute appreciably to 
the strangeness form factors via its coupling to the <fi. We return to this possibility in Section 
| (see also Ref. 



3 KK Intermediate State and Unitarity 

In order to determine KK contribution to the spectral functions, we need the matrix elements 
(N(p)\J N \K(k)K(k)) v(p) and (K(k)K(k)\ Jf ,\0). By expanding the KK -> NN amplitude in 
partial waves, we are able to impose the constraints of unitarity in a straightforward way. In 



doing so, we follow the helicity amplitude formalism of Jacob and Wick [22]. With A and A 
being the nucleon and antinucleon helicities, we write the corresponding S'-matrix element as 



(N(p,\)N(p,\)\S\K(k)KCk)) = 

i{27i) A 5\p + p-k-k){2n) 2 



6At 



t - Am\ 



1/2 



A,A|S(P)|00> 



where t = P 2 = (p + p) 2 and m# is the kaon mass. The matrix element (9, 0, A, A|S'(P)|00) is 
then expanded in partial waves as JIB], 



A, AI^C^IOO) = ^ f'^±l) //V A v! 



Ait 



o/xw? Q-i -(, 



(12) 



where T>^ u ,(a, f3, 7) is a Wigner rotation matrix with \i = A — A. The bj X define partial waves 
of angular momentum J. Because of the quantum numbers of the isoscalar EM and strange 
vector currents, only the J = 1 partial waves contribute to the spectral functions. Moreover, 
because of parity invariance only two of the four partial waves are independent. We choose 



6^ 2,1//2 and b^ 2 ' which fulfill the threshold relation [18 



1/2, -1/2 



(t) 



V2 b\ /2 ' 1/2 (t) 



t=im 2 N 



t=4m 2 N 



Using the above definitions, the unitarity of the ^-matrix, S^S = 1, requires that 

\b x / x (t)\<l, 



(13) 



(14) 



for t > Am 2 N . Consequently, unitarity gives model-independent bounds on the contribution of 

the physical region (t > Ami) to the imaginary part. In the unphysical region 

Am 2 N ), however, the partial waves are not bounded by unitarity. Therefore, we must rely upon 



7 



an analytic continuation of KN scattering amplitudes. This procedure is discussed in the next 
section. 

The second matrix element appearing in Eq. (0), (K(k)K(k)\j^\0), is parametrized by 
the kaon vector current form factor F": 



(0\j^\K(k)K(k)) = {k-k)»F«{t) 



(15) 



where a denotes EM or s and -F^(O) gives the corresponding charge (e.g., F^(0) = —1). 

From Eq. (|10D , the spectral functions are related to the partial waves and the kaon strangeness 
form factor by, 



ImF[ a) (t) = Re 



ImF^ a \t) = Re 



m N q t \ 


E 


4p? ) 


_\[2m N 


m N q t \ 


fe l/2,l/2 


*p 2 t J 





1/2, -1/2 



(*) 



m 



V2E 



b\ /2 > l/2 {t) 
(t) 



N ,1/2.-1/2 



where 



Pt 



E 



F a K {ty\ 

t/2. 



't/A-ml, q t = yft/l - m\ , and 
The corresponding spectral functions for the Sachs form factors follow from Eq. 

ImG^it) = Re 
ImG$(t) = Re 



Qt 



2V2t 



1/2, -1/2 



(t)F a K (ty 



(16) 
(17) 

(18) 

(19) 
(20) 



On one hand, Eqs. ( |TE| - PUD may be used to determine the spectral functions from experimen- 
tal data. On the other hand, one can impose bounds on the imaginary parts in the physical 
region by using Eq. ([14]) . Eqs. ( p^ - |20|) involve expressions of the type 

Re {b\ /2 ' ±1/2 (F^*} = \b 1 1 /2 ' ±1/2 \\F^\cos(6 1 -5 K ) = \b{ /2 ' ±1/2 \\F«\(l + lK ) , (21) 

where the phase correction ^ K is defined by Jk = cos(<5i — 5k) — 1, with 5\ and 5 k the 
complex phases of the b\' x and form factor, respectively. The experimental information on 7^ 
is incomplete. Since |1 + jk\ < 1, however, we can take 7^ = to obtain an upper bound on 
the spectral functions. In order to obtain finite bounds for the Dirac and_ Pauli form factors 
at the NN threshold, we build in the correct threshold relation for the 6^' A , Eq. fll3"|). This is 
necessary to cancel the 1/pf factor in Eqs. ( |TE| , |TTD . Strictly speaking, the relation holds only 
for t = Am 2 N . For simplicity, however, we assume this relation to be valid for all momentum 
transfers, as e.g. holds in the tree approximation of perturbation theory. Consequently, we 
have 



lmF[ a \t)\ < 
ImF^\t)\ < 



qt 



2\/2(v / t + 2m, 
m N q t 



2t(y/i + 2m t 



\m)\- 



(22) 



(23) 



The unitarity bounds for the Sachs form factors are obtained more straightforwardly by simply 
setting |6^' /2 ' ±1 ^ 2 



(*)|<1 in Eqs. QE9| 
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3.1 Kaon Form Factors 



The timelike EM kaon form factor, F^ M 



(t), has been determined from e + e~ — > KK cross 
sections. A striking feature of F^ M observed in these studies is the pronounced peak for t ~ m| 
|23fl . At higher values of t, oscillations at a much smaller scale are observed. A variety of 
analyses of F® M in this region have been performed |23| , |2"3|, and it is found that F^ M is 
well-described by a VMD parametrization: 



m 



1/ 



my — t — im v Tyf v (t) 



(24) 



where the sum is over vector mesons of mass my and width Ty and where fv{t) is some 
specified function of t. We use fv{t) = t/m v |[24|| . In nearly all analyses, one finds for the 
residues: C p ~ 1/2, C w ~ 1/6, and ~ 1/3. These residues may alternately be described in 
terms of the V'j and VKK couplings: Cy = gyxR/ ' fv , where f p ~ 5.1, f u ~ 17, and ~ 13. 
The strong couplings can be determined from r(0 — > KK) and SU(3) relations. 



To obtain F^(t) we follow Refs. JT5, 11], [2lJ and draw upon the known flavor content of the 
vector mesons. The p does not contribute to isoscalar form factors. To the extent that (i) the 
to and satisfy ideal mixing (|0) = —\ss) and |u;) = \uu + dd) / \f2) and (ii) the valence quarks 
determine the low-t behavior of the matrix elements (0|J M |V) one expects (0|s7 M s|u;) = and 
(O|s7 M s|0) = —3(017^^^10). It is straightforward to account for deviations from ideal mixing 



|4 nn 25: 



-Vq 



-V6 



sin e 



sin(e + On 



cose 



cos(e + 6*0 



-0.2 



(25) 



where the s superscript denotes the residue for the strangeness form factor, #0 is the "magic" 
octet-singlet mixing angle giving rise to pure uu + dd and ss states and e deviations from ideal 
mixing. From Eqs. ([24]) and (^) we observe that the time-like kaon strangeness form factor is 
dominated by the 0(1020) resonance. We note that the flavor rotation of Eq. (^) only gives 
the relative size of the u and contributions but does not lead to the correct normalization for 
at t = which must be enforced by hand. 

In Fig. [I] we plot F^ as given by Eqs. Q2"4] ) and (^) and compare it with a simple VMD model 
|18| and the Gounaris-Sakurai (GS) parametrization for F n [26] with all p parameters replaced 



by the corresponding ones for the 0. We observe that the GS and VMD forms reproduce the 
essential features of F* as determined from e + e~ data and standard flavor rotation arguments. 
Since F^ is needed for t > ^m 2 K1 the (^-contribution which gives rise to the bump around t m 
0.6 GeV 2 in Fig. |I] is negligible. In comparison to the strong peak, the small scale oscillations 
at higher-t have a negligible impact as well. When computing the leading strangeness moments, 
we find less a than 10% variation in the results when any of these different parametrizations 
for Fx is used. In short, any model parametrization of showing the peak at the mass and 
having the correct normalization, F£(0) = —1, may be used for the purpose of studying nucleon 
strangeness. The pointlike approximation, F^ = —1, however, misses important resonance 
physics. Throughout the remainder of this work, we will use the GS parametrization. 
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10 

0.0 0.4 0.8 1.2 1.6 2.0 

t [GeV 2 ] 

Figure 1: Different model parametrizations for F^. Full line shows simple VMD model with the 
only, dashed line shows the flavor rotated VMD model including the u and <f), and dash-dotted 
line shows GS parametrization. Dotted vertical line indicates KK threshold. 

4 Analytic Continuation 

To obtain the b\' for Am 2 K < t < Am 2 N we analytically continue physical amplitudes into the 
unphysical regime. The analytic continuation (AC) of a finite set of experimental amplitudes 
with non-zero error is fraught with potential ambiguities. Indeed, AC in this case is inherently 
unstable, and analyticity alone has no predictive power. Additional information must be used 
in order to stabilize the problem, as we discuss below |27j . 

In order to illustrate these issues and the methods we adopt to resolve them, we first briefly 
review the kinematics of KN scattering (cf. Ref. P§1). It is useful to consider the s-, u-, and 
t-channel reactions simultaneously,^] 

(a) s-channel: K( qi ) + N( Pi ) -> K(q f ) + N(p f ) , (26) 

(b) u-channel: K(-q f ) + N( Pi ) -> K(- qi ) + N(p f ) , 

(c) t-channel: K(-q f ) + K( qi ) -> N(- Pi ) + N(p f ) , 

where the four-momenta of the particles are given in parentheses. In this notation the crossing 
relations between the different channels are immediately transparent. The three processes can 
be described in terms of the usual Mandelstam variables 

s = (Pi + qi) 2 , u = (pi- q f ) 2 , and t = (<& - q f ) 2 . (27) 

4 We do not consider K°N scattering data in this analysis. In the following, we write K and K for K + and 
K~, respectively. 
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The invariant matrix element for the KN scattering process has the structure, 



M = u(pf) 



(2* 



and the isospin decomposition of the invariant amplitudes A and B reads 



A(s,t) = A + (s,t)+A~(s,t)(T N -f K ) } 
B(s,t) = B + ( Sl t) + B-(s,t)(T N -T K ) . 



(29) 



We also need the A and S pole contributions to A^ and B^' which are given by 



^pole( S ' U ) 



E9kny I \ 
— — {m Y -m N ) 



1 ±_L 



y=A,s 



^ u — m\ s — m Y J ' 



E 

y=A,s 



9kny 



u — m Y s — m Y t 



(30) 




Figure 2: Mandelstam plane for KN scattering. Physical regions are marked by dashed areas. 

It is instructive to display the range of s, u, and t in the Mandelstam plane, shown in Fig 
0. The physical regions of the three reactions do not overlap, and the invariant amplitudes 
simultaneously describe all three processes. The physical values of the invariant amplitudes are 
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obtained when the Mandelstam variables are taken in the corresponding ranges. In order to 
carry out the dispersion integrals of Eqs. ([?J |]), we require the 6^' A along the t-channel cut, 
indicated by the gray shaded area in Fig. ||. 

We begin with experimental K N amplitudes in the s-channel region and employ the method 
of backward DR to obtain the unphysical amplitudes along the t-channel cut. The backward DR 
method has has been used successfully for a similar continuation of irN scattering amplitudes 
|29| and as a consistency test for different KN phase shift solutions in the 1970's [0. Only 



recently has a continuation of the KN amplitudes become possible due to improvements in the 
data base over the last three decades. A continuation of KN scattering amplitudes has also 
been performed on the basis of hyperbolic DR . In the latter analysis, however, the t-channel 
helicity amplitudes have been parametrized by sharp resonance poles, an a priori assumption 
we seek to avoid in the present analysis. We use the recent KN phase shift analysis of the 



VPI-group |32j as experimental input. 



4.1 Overall Strategy and Problems 

Although an AC is uniquely defined from a continuum of points, this is not the case for a finite 
set of points lying within a certain error corridor. Consequently, the procedure for obtaining 
amplitudes outside the range of the given data is unstable. In practice the problem is stabilized 
by restricting the admissible solutions. This is achieved by using a priori information on the 
solution from experiment or theory as "stabilizing levers" p7f . In the remainder of this sub- 



section, we give an overview of our general procedure for making this stabilized continuation 
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The first important observation is that in the backward direction, the invariant amplitudes 
in the s- and t-channel coincide, i.e. f s (t,9 s — ir) — ft{t,6t = n) = f(t). Here / is a 
generic invariant amplitude and 9 S and 9t are the CM-scattering angles in the s- and t-channel, 
respectively. Thus, the continued amplitude f s (t,9 s = 7r) in the unphysical regime can be 
identified with the unphysical t-channel amplitude, ft(t,9 t = ir). 

Second, we work on a Riemann sheet where only the singularities from the s- and t-channel 
reactions are present. The backward amplitude / then possesses a left-hand cut from zero to 
minus infinity stemming from the s-channel reaction. The lowest right-hand cut, which is due 
to the 3tt intermediate state in the t-channel, runs from 9m 2 to infinity. The two cuts do not 
overlap and we can write an unsubtracted DR for /. The amplitude is known from experiment 
on the left-hand cut in the region t p < t < and we need its values on the right-hand cut where 
/ is related to the KK — > NN partial waves. In particular, we are interested in the unphysical 
region 4m 2 K < t < 4m^ where the partial waves are not bounded by unitarity. 

Following our strategy to include as much a priori information as possible, we do not continue 
the amplitude / itself but rather a related function, Af, the so-called discrepancy function. In 
Af the analytically determined pole terms for /, as well as the experimentally known behavior 
of the amplitude, have been subtracted out. This subtraction is carried out in such a way 
as to remove the portion of the left-hand cut lying in the range t p < t < where the phase 
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shift analyses are available. Thus, the use of Af presents two advantages. First, the original 
problem which entailed an AC from the boundary of the region of analyticity to another point 
on the boundary (B — > B) has been transformed into one in which the continuation occurs 
from the interior of the analyticity domain to the boundary (/ — > B). Second, the well-known 
pole terms can be continued explicitly, and only non-pole parts which remain in Af must be 
continued with other methods. In the end, / is easily reconstructed from Af. 

The continuation of Af is carried out by means of a power series expansion. Such an 
expansion, however, converges only up to the nearest singularity. In our case this point lies at 
the three-pion cut, well below the region of interest. We circumvent this problem by using a 
conformal mapping of the complex t-plane. First we symmetrize the cuts to lie along (— oo, —R] 
and [R, oo). We then map the cuts for Af onto the circumference of the unit circle in the w- 
plane and the rest of the t-plane into its interior using 

. . 2Rw , . 

t(w) = — — . 31 
1 + «r 

Now it is possible to use a Legendre series in w to perform the continuation. In principle, the 
series can be fit to Af in the region t p < t < and the fit then be evaluated in the region of 
interest. In practice, we fit to the function 

Ag(w) = \t\- a Af(t(w)) (32) 

instead of Af. In so doing, we explicitly enforce the high-energy behavior of the amplitude /: 
f(t) ~ \t\ a for t — > oo. The parameter a may be obtained from Regge analyses. In addition, 
the explicit inclusion of the high-energy behavior stabilizes the continuation by suppressing os- 
cillations in the solution. In fact, a can also be thought of as an auxiliary parameter introduced 
just for this purpose. The final results depend only weakly on the precise value of a. In fact, we 
evaluate the stability of the continuation by studying the variation of the results with a as well 
as with the number of terms included in the Legendre series. The technical details of the AC 



procedure are given in Ref. l28fl . In the next subsection, we directly start with the continued 
backward amplitudes. 



4.2 Results 

First we introduce the invariant amplitude which in the backward limit is related to 
and flW via 

jX±)(t) = + m N J^ t B^\t) . (33) 

In order to extract the 6^' A , we need B^ and in the t-channel unphysical region. In 
Fig. |3|, we show the results of the analytic continuation for FW = _ F^ e and B^ = 
B^ — -Bpo/ e up to t = 14 vr? K . The continued amplitudes vanish as t — » oo. We do not 
expect the power series expansion to be credible in a larger region than the one in which 
its coefficients have been determined. Since the experimental amplitudes are given in the 
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Figure 3: Analytic continuation of the invariant backward amplitudes = — F^ e (with 



n 



6) and B(+) = B« - B^ e (with n = 5) for a = -1.2. 



interval — 9>m K < t < 0, we trust our continuation only up to t ~ %m K which covers half of 
the unphysical region. Fortunately, for purposes of analyzing nucleon form factors, the KK 
threshold region dominates. The kaon form factor F£(t), which multiplies the b x,x in the spectral 
functions amplifies the low-t contribution and suppresses that from the region 8m 2 K < t < 4m 2 N . 
We plot only the results for the amplitudes without their pole part. Since the pole part is 
exactly known, it can be continued explicitly and is added to the remainder obtained from the 
AC in the end. 

Ultimately, we require the J = 1 partial wave projections of the backward amplitudes which 
can be expanded in partial waves as 



F (+ \t) 



57T 



E 



(t)P'A 



J + 1 ft 1 /-' - ^MP' 

Qt ~! y/j(J+l) 

AnVi ^ 1 



(34) 



Vtlt T-n * 



j=o 



Because the AC gives the invariant amplitudes only in the backward direction, it is non-trivial 
to extract the partial waves. The expansion can be carried out separately for / and / po i e - This 
is useful because the sum in Eq. (^) converges much faster for the non-pole part than for the 
pole part f29fl. Consequently, we exploit the faster convergence for the non-pole part and add 
the exactly known pole term projections after the bi' x have been isolated. 

Since the amplitude is known only for one value of 6 t , a separation of the 5^ 2,±1 ^ 2 from the 
bj'^ relies on several additional observations. First, each \bj ,x \ — > 1 or smaller as t — > Am 2 N 
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because of unitarity. The only significant deviation from this trend for 4m 2 K <t< Am 2 N occurs 
via resonant enhancements of the partial waves. The lightest J > 1, I = resonance having a 
non-negligible branching ratio to the KK state is the /2(1270), whose mass lies near the upper 
end of the range of validity of the AC. Moreover, in Ref. |31] no evidence for J > 2 resonance 
effects close to the KK threshold was found. Consequently, we truncate the expansions in Eq. 
(531) at J = 1. In order to test the validity of this truncation, we also perform an analysis with 
a model for the resonant J = 2 partial waves included. As discussed below, our results are 
essentially unaffected by this inclusion. Since the width of the /2(1270) is ~ 185 MeV, we may 
expect some small, residual contamination of the J < 1 partial waves due to our truncation. 
Fortunately, the presence of F£(t) in Eqs. (fTT)| - |2T)|) protects the spectral functions. The kaon 
form factors are peaked in the vicinity of the 0(1020) and strongly suppress contributions from 
t > ml 



The remaining separation between the S- and P-waves in 

£(+) 

may be performed by drawing 

on the work of Refs. [ 30| and |H] , where backward and hyperbolic DR have been used to analyze 
the KK — ► NN amplitudes under the assumption that the helicity amplitudes are dominated 
by sharp resonance poles. In the S-wave, one finds a resonance close to the KK threshold 
having a 22% branching ratio to KK, namely the /o(980). Therefore, we use the resonant 
6q //2 ' 1/ ' 2 amplitude from Refs. |30], [TT[] and subtract it from our result.^ Both analyses use the 
following form for the amplitude 



(1/2,1/2 



(t) 



2q t 



\/ip t m 



t 



(35) 



and obtain approximately the same result for the residue r+ . We use r+ 
[3~T| and extract the b{ from the continued invariant amplitudes using 



;-24.10±6.65) m 



1/2, -1/2 



(t) 



- 1/2 (t) 



12tt 



pole 



(36) 



b{ /2 ' 1/2 (t) 



i(^ (+) (*) + 6S Al/2 ) + ^ 1/a (*) 



pole 



with p. 



t/A. Note that the A and S pole term projections for even values of J vanish 



for the isoscalar amplitudes in the t-channel. As a consequence, the pole term projection for 
J = does not appear in Eqs. (|36|). The J = 1 pole term projections are 



b\ /2 >- 1/2 (t) 
b{ /2 > 1/2 (t) 



pole 



pole 



Y=A,E 



(37) 



E 

y=A,s 



9kny 



m N q t 



+ 



Sip. 



(m Y - m N ) Qi{-i£ Y ) 
(2Q 2 (-z£ y )+QoH£0) 



5 We can not use the results of Refs. [^0|, |3lJ directly because the 6^' A were assumed to be dominated by a 
single effective u> pole and the explicit effect of the <f) was a priori excluded. 
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where the Qi, i = 0,1,2 are Legendre functions of the second kind and £ y is given by 

t - 2m 2 K + 2{m\ - m 2 N ) 



(38) 



The introduction of a width in Eq. (^) changes the extracted amplitudes b\j 2, ±1//2 only slightly 
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Figure 4: Comparison of the results for the fr}^ 2 ' 1 ^ 2 (upper row) and &| /z '~ i/z (lower row) for 
different order of the Legendre series n and a = —1.2. Solid lines show |, dashed 

lines show i?e&J //2 ' ±1//2 , and dotted lines show Im6^ 2,±1 ^ 2 .The continuation is expected to be 
trustworthy up to t w 2 GeV 2 . 

and has almost no impact on the application to nucleon form factors. The same observation 
applies for the case when the J = 2 partial waves from Ref. [31] are subtracted as well. 
Therefore, we discard the J > 2 partial waves and proceed accordingly. 

In the following, we discuss the sensitivity of the amplitudes to the order n of the Legendre 
series and the high-energy parameter a. In Fig. || we display the full &J' /2 ' ±1 ^ 2 for different 
values of n. We show their real and imaginary parts as well as their absolute value. It is clearly 
seen that bY 2 ' 1//2 is almost independent of n, whereas b\^ 2 ' ~ 1 ^ 2 shows some variation. From the 
quality of the fits to the discrepancy functions, we choose n = 6 for F( + \ whereas we take 
n = 5 for fiW. This choice corresponds to taking the minimum n which gives a satisfactory 



1/2, -1/2 



fit in order to minimize the amplification of experimental noise |27| , p5| . The dependence of 
our results on the asymptotic parameter a is similar. Although a is a physical parameter, its 
determination is model dependent. We have varied a from —0.2 to —5.0 to test the sensitivity 
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1/2 1 /2 

of the AC to a. Similar to the dependence on n, the b x ' are almost independent of a, 
whereas the fo^ 2 '" 1 ^ 2 show some variation. For a — > our analysis becomes unstable as we 
expect since arbitrary oscillations in the AC are no longer suppressed. Furthermore, for a > 
the DR does not converge and the corresponding results are meaningless. For the final results 
we take a = —1.2 from the Regge- model fit of Ref. [Q. From the dependence on n and a, we 



expect our results for b 1 ^ 2 ' l ^ 2 to be more reliable. This conclusion is supported by the fact that 
\F\ 2 is directly related to the backward cross section. Hence, F is particularly well determined 
by the data and errors in the phase shift analysis cancel in the reconstruction of F |29|, |30fl . 



As shown in Fig. £|, a clear resonance structure at threshold is seen in b\ 2 , which 
presumably is the <\> resonance. This resonance is not observed in 6^ 2, 1 ^ 2 . There are two 
possible explanations for this fact: 

1. The AC for B and consequently b^ 2 ' 1 ^ 2 is not sufficiently well determined for the reasons 
explained above. 

2. The vector and tensor couplings of the <fr meson to the nucleon are approximately equal and 
have opposite signs. Parametrizing the resonant part of the J = 1 partial waves by the 0, we 



have 31 



, 2^2 G^K , y T \ 

L - — o a \ (f>NN ' ^(pNN J 1 \ OJ ) 



An 

2m N G^ KK ( ' v m n T 

" ' rhNN + . o " '^ 



+ - 3 An \ "* NN ^ Am? N * NN ) ' 
where the T± are the residues corresponding to 5^ 2,±1//2 ( c f. Eq. (^)). GX NI ^ = —G^ N ^ then 



' cj>N N " <j>NN 

leads to a resonance in b^ 2 ' 1 ^ 2 but not in b\ 2 . 

Following scenario 2 above, we determine GX N ^ and G^ N ^ from the resonance structure in 
^/ 2il / 2 _ ^y e ^ e region around this structure using the J = 1 analog of Eq. fl35|) , but 
including a finite width, r^. We obtain G^ N ^ = 7.4 using the Particle Data Group value 



for |36| and GY N ^ = 9.6 when is allowed to be a fit parameter. These values are 



comparable with the value GX N <j = 9.2 obtained from the VMD parametrization of F[' 0) [[14 



' 4>NN 

Our value for G^ N ^ = —GY N ^, however, has a larger magnitude than obtained in Ref. JT1 



The possible reasons for this difference are discussed elsewhere ||37|| . We emphasize, however 



that our determination of the <pNN couplings relies only on the observed resonance structure in 

1/2 1/2 

the strong amplitude b-y ' and not on a VMD ansatz for the isoscalar EM form factors. We 
also note that the interpretation of the by ,x in terms of (pNN couplings is inconsequential for 
the computation of the F i , since the amplitudes themselves - rather than a parametrization 
of them - are used in the dispersion integrals. 
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5 KK contribution to Nucleon Form Factors 

We now use the b±' x to evaluate the contribution of the unphysical region to the dispersion in- 
tegral for the strange and isoscalar nucleon form factors. We consider the electric and magnetic 
Sachs radii and the anomalous magnetic moment. Drawing upon Eqs. (0, H), we write down 
subtracted DR's for the Sachs radii, 

4m 2 r°° ImG a (t') 
P-=~^ / 2 dt' 1 -^^!, (40) 



K 



where i = E,M and a denotes the isoscalar EM or strangeness channel. Before evaluating Eq. 
(f40|), we discuss the qualitative features of the spectral functions for p E , p a M) and n a . Since for 



our purpose the difference between the EM and strange kaon form factors is essentially given 
by the normalization, the qualitative features of the strange and EM spectral functions are the 
same. Consequently, we limit the following discussion to the strange spectral functions. 

The KK content of the spectral function for p s E is displayed in Fig. [5]. In Fig. |5|a, we 
compare two scenarios for b^ 2 ' 1 ^ 2 : (I) the Born approximation in the nonlinear cr-model (BA) 
and (II) the AC from the previous section. For both scenarios a pointlike kaon form factor 




Figure 5: KK contributions to the spectral function for p s E . Short dashed curve (a) gives 0(g 2 ) 
result of scenario (I) discussed in the text. Solid curve [(a) and (b)] gives result for scenario (II). 
Long-dashed curve gives unitarity bound for t > 4m^ using a pointlike (scenario (II)) and 
GS parametrized F£ (scenario (III)). Dashed-dotted curve (b) gives all-orders spectral function 
(scenario (III)). Dotted vertical line indicates physical NN production threshold. 

F£(t) = 1 has been used. Although the pointlike form factor is unrealistic, using it in this 
context allows to illustrate separately the effects of form factor and scattering amplitudes. 
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We include scenario (I) because of its correspondence with a number of model calculations as 
well as CHPT. The scenario (I) spectral function contains only contributions to 0(g 2 ), where g 
is the scale of the strong hadronic couplings. It constitutes the DR form of a one-loop calculation 
containing a kaon and strange baryon intermediate state and a current insertion on the kaon line 
|18| , p0| . A variety of model calculations (see, e.g., Refs. ||, ||) have been performed under the 
assumption that such amplitudes dominate the strangeness form factors and that a truncation 
at 0(g 2 ) gives a reliable estimate of the scale and sign of the kaon contribution. In ChPT 
analyses of both the strange and isoscalar EM moments, the non-analytic (in quark masses) 
parts of the same amplitudes are retained and added to the appropriate low-energy constants. 
In the case of p s E , for example, the leading-order non-analytic contribution is singular in the 



chiral limit 11 



P 



non-anal 




In 4 



(41) 



where F and D are the usual SU(3) reduced matrix elements, A x = 4:7if n gives the scale of 
chiral symmetry breaking, and p is a renormalization scale. The chiral singularity of Eq. (fTI|) 
arises from the kaon-nucleon one loop graph with the strange vector current inserted on the 
kaon line. The equivalence of this loop calculation with the scenario (I) for the DR analysis 



|18| , |20| , [211 guarantees that the DR computation includes p non - ana \- 

At leading order ChPT, both p s E and k s also contain low-energy constants which cannot be 
obtained from existing data using chiral symmetry JTTJ . For the magnetic radius, however, only 
the non-analytic one-loop terms contribute at 0(p 3 ) |L2 |. In principle, the low-energy constants 
include higher-order effects (in g) not contained in the non-analytic loop contributions. A 
comparison of the scenario (I) spectral function with the full spectral function - containing 
contributions to all orders in g - allows us to determine the importance of these low-energy 
constants and more generally to evaluate the credibility of one-loop predictions. 

For both scenarios (I) and (II), we show the upper limit on the spectral function generated 
by the unitarity bound on b 1 / 2 ' 1 ^ 2 . As observed previously in Ref. |18| , and as illustrated in 
Fig. [5]a, the BA omits these rescattering corrections and consequently violates the unitarity 
bound by a factor of four or more even at the NN threshold. This feature alone casts a 
shadow on predictions which rely solely on one-loop amplitudes. In principle, the low-energy 
constants of ChPT correct for the unitarity violation implicit in the one-loop contributions. The 
counterterm-free 0(p 3 ) prediction for p s M , however, does not include this unitarity correction. 

The curve obtained in scenario (II) indicates the presence of a peak in the vicinity of the 
0(1020) meson, reflecting the presence of a KK <-» <fi resonance in b 1 ^ 2,1 ^ 2 . This structure 
enhances the spectral function over the 0(g 2 ) result near the KK threshold. As t increases from 
4:m 2 K , the spectral function obtained in scenario (II) falls below that of scenario (I), presumably 
due to KK rescattering which must eventually bring the spectral function below the unitarity 
bound for t > Am 2 N . The fact that the AC is approximately constant above t > 2 GeV 2 and 
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also violates the unitarity bound, indicates that the AC and partial wave separation cannot be 
trusted close to the NN threshold. 

In Fig. [5|b we plot the spectral function for a third scenario (III), obtained by includ- 
ing a realistic (empirical) kaon form factor. We use the GS parametrization, although other 
parametrizations for F^, such as a simple (^-dominance form, yield similar results as has been 
discussed above and elsewhere |Tj|. The corresponding spectral function is shown in Fig. ||b, 
and compared with the spectral function for a pointlike F£. The use of the GS parametrization 
significantly enhances the spectral function near the beginning of the KK cut as compared 
with the pointlike case, while it suppresses the spectral function for t > 2 GeV 2 . Consequently, 
the full spectral function is dominated by the low-t region where the AC for the b x ' x is reliable. 
The t > 2 GeV 2 region gives a negligible contribution, even though the b x ' x are too large in this 
range (see Fig. |5|b). The error associated with this region is correspondingly negligible. We 
emphasize that the one-loop model predictions miss entirely the resonance enhancement of the 
KK contribution. 

The corresponding scenarios for the magnetic Sachs radius p s M are shown in Fig. [| In 
contrast to 6^ 2,1 ^ 2 , 6^ 2 ' -1 ^ 2 shows no resonance structure at all. The possible explanations for 
this behavior have been discussed in the previous section. Here we take our result at face value 
and develop the pertinent consequences for p a M . In Fig. ^|a we compare the results of scenario 




Figure 6: KK contribution to spectral function for p s M . Curves are as in Fig. |5j. 

(I) BA and (II) AC for a pointlike coupling of the kaon to s7 M s. The AC spectral function stays 
clearly below the BA result for all t, presumably due to rescattering effects included in the AC. 
Similar to the electric case the unitarity bound is also violated in scenario (II), although the 
violation is about a factor two weaker than in the BA. Again we conclude that the continued 
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b^ 2 ' 1,/2 is not trustworthy close to t — 4m 2 N . The suppression of the spectral function for 
t > 2 GeV 2 by the GS form factor in scenario (III) is displayed in Fig. |(|b. 

1/2 1/2 

The absence of the resonance structure in 6/ ' suggests that the integrated spectral 
function for scenario (III) gives an upper bound, rather than a firm prediction, for p a M . The 
reason has to do with the relative phases of b± ' and F£(t). Consider the strangeness case 
with a simple VMD parametrization of F^(t) (other parametrizations are similar): 

2 

F s K (t)= 2 , *. r ■ (42) 

As input to Eqs. (28-32), we require 

Rebf X F°* = (Re b x ' x )(Re F S K *) + [Im b x ~ x )(Im F S K *) , (43) 

where Re b x,x and Im b x,x are shown in Fig. As t crosses the resonance in the vicinity of 
rr&, Re F^* changes sign while ImF^* does not. Instead, the latter reaches its peak value of 
magnitude ~ m^/T^. From Fig. |], we observe that neither Reb x,x nor Imb x,x undergoes a 
phase change around t = rr&. Hence, when integrated across the resonance, the contributions 
to the integral from (Reb x ' x )(ReF^*) change sign, leading to substantial cancellations. The 
contributions from (Imb X ' X )(ImF^*), on the other hand, do not change sign, and no cancel- 
lations occur. In the case of bY 2 ' 1 ^ 2 , it is Imfe^ 2 '^ 2 which contains the resonance structure. 
Consequently, the integral for p s E is dominated by the resonating imaginary parts of b^ 2 ' 1//2 and 
F^, which remain in phase. The continuum contribution is suppressed by the relative phase 
change of the real parts. 

The situation for p s M is different because b^ 2, 1 ^ 2 displays no resonant behavior. Again, 

i/o ]_ /2 

the contribution to the integral from [Rebi ' )(ReF^*) is suppressed by cancellations from 
the phase change across the resonance. |(Im6^ 2 ' 1 ^ 2 )(ImF^*)\ <C |&i^ 2, 1 ^ 2 F^*\, however, and 
its precise magnitude is sensitive to the parameters a and n. Consequently, the integral of 
Re(b\ 1 F^*) is rather uncertain. We are confident, therefore, in quoting only an upper 
bound for p s M , obtained by integrating \b± ' F^ \, which varies only gently with a and n, 
and using ^ K = 1 (Eq. (|2T|)). As we note below, even this upper bound is nearly twice as small 
as the result obtained from scenario (I). 

To determine the anomalous magnetic moment K, a , we require an unsubtracted DR and turn 
to Fg . Eq.([7p reduces in the limit t = to 

K a = F^) = - f°° dt' lmF2itl) . (44) 



71 Jim 2 K 



with IinF^ given in Eq. (|T7|). However, an additional comment is in order. ImF£ depends on 
both b x,x rather than on one as for the radii. In order to guarantee a finite spectral function at 
the NN threshold, the two b x,x must fulfill the threshold relation b^ 2 ' 1 ^ 2 = \f 7 LbY 2 ' ^ 2 . Our 
AC, however, is not reliable at t = 4m^ and does not obey this relation. Therefore, we replace 
ImF£(t) by ImF$(t = 8m 2 K ) for 8m 2 K < t < Am 2 N . Doing so leads to an upper bound for the 
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Scenario 


Moment 


A 9 ^ , ^ A 9 

Am K <t< Am z N 


A 9 ^ , 

Am z N < t 


Total 






0.25 


0.10 


0.35 


0(g 2 ) 


Pm 


0.30 


0.05 


0.35 




K S 


-0.07 


-0.07 


-0.14 




Pe 


0.98 


0.01 


0.99 


AC/GS 


Pm 


0.17 


0.01 


0.18 






-0.41 


-0.01 


-0.42 



Table 2: Kaon cloud contribution for p s E , p s M , and k s for the scenarios discussed in the text. 
Second and third columns give contributions to the dispersion integral from theunphysical and 
physical region, respectively. In both scenarios, the unitarity bound on the b\' x is imposed for 
t > 4m 2 N . To convert p s to (r 2 ) s , multiply p s by -0.066 fm 2 . 

spectral function. Essentially the same qualitative observations as for p E apply to n a because 
the spectral function is dominated by the resonance in b 1 ^ 2, 1 ^ 2 in both cases. Specifically, the 
contribution from t > 8m 2 K is negligible. 

The numerical consequences of our analysis are indicated in Table |2], where we give results for 
the leading strangeness moments. The first three lines give the 0(g 2 ) "kaon cloud" prediction 
for p s E , p s M , and k s (scenario I). The last three lines give the results when the GS form factor 

1/2 | "I /2 

and analytically continued fit to KN amplitudes for b x ' are used (AC/GS). In all cases 
the unitarity bound on b^ 2 ^ 1 ! 2 is imposed for t > Am 2 N . The AC/GS results were obtained 
by extending 6^ 2,1 ^ 2 to 4m 2 . Although the continuation can be trusted only for t<8m 2 K , and 
although the continued amplitudes exceed the unitarity bound for t — *> Am 2 N , the GS form factor 
suppresses contributions for t > 8m^, rendering the overall contribution from 8m 2 K < t < 4m 2 N 
negligible. The overall sign of the product of F^* and the }^ 2,±x l 2 has been determined by 
fits to EM form factor data f3~?| . However, the sign of p s E is well determined from the phase of 
b^ 2 ' 1 ^ 2 and F^* at the 4> peak. 

From Table ^|, we observe that the use of a realistic KK spectral function increases the kaon 
cloud contribution to p s E by roughly a factor of three as compared to the 0(g 2 ) calculation. 
Moreover, the result of the all-orders computation approaches the scale at which the proposed 
PV electron scattering experiments are sensitive, whereas the 0(g 2 ) calculation (e.g., one- 
loop) gives a result which is too small to be seen. This observation depends critically on the 
presence of the resonance structure near t = m 2 . in both b 1 ^ 2 ' 1 ^ 2 and F^ (Fig. |). Its absence 
from either the scattering amplitude or kaon form factor would lead to a significantly smaller 
magnitude for p E . 

Because b 1 ^ 2, 1 ^ 2 does not display any resonance structure, no enhancement in p s M is ob- 
served. In fact the 0(g 2 ) result is a factor of two larger than in the AC/GS case. We empha- 
size that the AC/GS result for p s M constitutes an upper bound, given the phase uncertainty 
discussed above. It is noteworthy that this bound \p s M \ = 0.18 is not consistent with the range 
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given by ChPT to order 0(p 3 ), p s M = 2.44 . . . 9.05 0. Although both the DR and ChPT cal- 
culations of p s M include only the KK contribution, a large discrepancy exists between the two 
methods. We suspect that the problem lies in the truncation at 0(p 3 ) in ChPT, for the follow- 
ing reasons. First, the 0(g 2 ) spectral function contains no rescattering corrections required for 



consistency with the unitarity bound. When integrated as in Eq. (|4*0|), this spectral function 
yields the value for p s M given in Table [|. Second, the 0(p 3 ) calculation in ChPT includes only 
the non-analytic contributions to this integral. The latter yield a result that is an order of 
magnitude larger than the result obtained for the full 0(g 2 ) integral. Given that this is already 
a factor of two larger than the upper bound obtained from the all orders DR result, we suspect 
that the 0(p 3 ) ChPT result omits crucial higher-order rescattering effects. Presumably, these 
effects are contained in the 0(p A ) Lagrangian, which contains a term of the form 

£(4) = -^^^ 2 F^, (45) 

x 

where A x = 4^/^. In effect, the AC/GS result of Table represents a determination of the 
kaon rescattering contributions to 6^' A . Unfortunately, ChPT cannot make counterterm free 



predictions for the other moments jrl]] an d the problem remains unresolved. 

In the case of k s , we observe an increase by a factor of three over the 0(g 2 ) approximation. 
However, due to the presence of both 6} //2 ' 1//2 and 6^ 2,_1//2 in ImF 2 and the phase-related 
uncertainty associated with 6 1 , we consider the all-orders AC/GS value for k s to be an 
upper bound. Given the size of the experimental errors, this bound is not incompatible with 
the SAMPLE results (see Eq. (@)). 

It is interesting to note that only the normalization of G S M at q 2 = but not its g 2 -dependence 
receives a resonance-enhanced KK contribution. This feature can be understood by observ- 
ing that the resonance-enhanced partial wave 6} //2 ' 1//2 enters ImF\ and Jmi 7 ^ with the same 
coefficient but opposite sign (cf. Eqs (§], [16], [L7]))- 

In Fig. 0, we plot the quantity extracted from forward angle parity-violating electron scatter- 
ing experiments on the proton, G s e + t(1 + k p )G s m , as a function of the momentum transfer t. At 
t = 0.48 GeV 2 , one can compare the result of the HAPPEX collaboration (see Eq. (0)) with the 
KK contribution: \G S E + 0.39G S M \ K ^- = 0.053. The two values agree within the experimental 
error bars.[] 

From a qualitative standpoint, the KK spectral content of the isoscalar EM and strangeness 
form factors are similar. The numerical significance, however, differs between the two cases. 
In Table BL we give the KK contribution to the leading isoscalar EM moments. We compare 
the KK contributions for the AC/GS scenario with the total "experimental" values from the 
dispersion theoretical analysis of Ref. |14| and the <ft pole contribution alone.f] Evidently, the 
g 2 -dependence of the isoscalar EM form factors is determined by states other than \KK). Its 
contribution to k^ i=0 \ however, is considerable. The VMD analyses of Refs. 116], O indicate 



a significant contribution to the leading moments from the 0. The results of the VMD treatment 
and our present study are compatible only if other intermediate states having t < vnh contain 



considerable (p strength. A preliminary consideration of this possibility is given in Ref. P5 ], 
where the role of the 6 resonance in the 3ir channel is discussed (see also Ref. |3 



6 Note that we quote only the absolute value because of the phase-related uncertainty in G S M . 
7 If a phenomenological EM form factor of the kaon (cf. Eq. ((2^)) is used the numbers in the second line of 



Table | are reduced by 20%. 
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Figure 7: KK contribution to G S E + r(l + K P )G S M as function of the momentum transfer £. The 
contributions from G% and G^. enter with opposite signs. 

6 Other Contributions and Conclusions 

The role of continuum and resonance effects in the isoscalar EM and strangeness form factors 
appears considerably more complicated than in the case of the isovector EM form factors 
whose spectral functions are dominated by a combination of p resonance and uncorrelated 
7T7T continuum. In the present study, we have continued our previous efforts |0| |25|, [TjJ to 
determine the corresponding picture for the isoscalar EM and strange vector current spectral 
functions. We have focused on the KK contribution for two reasons: (i) the availability of 
scattering data afford us with the least model-dependent determination of this contribution to 



Scenario 


(1=0) 

Pe 


(7=0) 

Pm 


„(7=0) 


Exp't 


-4.55 


-3.97 


-0.06 


AC/GS 


-0.50 


-0.09 


0.21 


pole |g 


2.21 


2.87 


0.15 



Table 3: Leading moments of the isoscalar EM nucleon form factors. First line gives "exper- 
imental" values from the dispersion analysis of Ref. [Q. Second line shows KK contribution 



in AC/GS scenario. Third line gives contribution of pole alone in the pole fits of Ref. [IE 
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all orders in the strong coupling, and (ii) the OZI rule has prompted a number of strange form 
factor calculations assuming this state to give the dominant contributions. We find that 

(a) the KK contribution to the isoscalar EM and strangeness electric spectral functions is 
significantly enhanced by the presence of a ^-resonance in 6^ 2,1 ^ 2 and F®. 

(b) there exists no evidence for such a resonance in the b^ 2 ' partial wave. 

(c) the resonance affects only the normalization of G a M but not its g 2 -dependence. 

(d) results (a) and (b) can be reconciled with a simple ^-resonance model of NN — > KK if 
the vector and tensor (f>NN couplings have roughly equal magnitudes and opposite signs. We 
obtain a value for G^ N ^ in agreement with the VMD analyses of the isoscalar EM form factors 
flUl . Our value for G^ N ^, however, is larger in magnitude. 

(e) the KK contribution to the magnetic radius p s M is significantly smaller than the value 
obtained at 0(p 3 ) in ChPT. 

(f) the KK contribution to the sub-leading g 2 -dependence of the isoscalar EM moments is 
small. 

The result (f) implies that consideration of other intermediate states is essential to a proper 
description of the isoscalar EM and strangeness spectral functions. In this respect, the calcula- 
tions of Refs. || [3^ are suggestive, indicating the possibility of cancellations between different 



contributions as successive higher-mass intermediate states are included. Moreover, as discussed 
in Refs. [^, |37| , contributions from light, multi-meson intermediate states may be just as large 



as that of the KK state. Nevertheless, our study of the KK state provides several insights 
into the the treatment of these contributions. In particular, we are able to understand the con- 
nection between continuum and resonance contributions to the isoscalar EM and strangeness 
form factors and to evaluate the credibility of other approaches used in computing them. In- 
deed, perturbative calculations which truncate at 0(g 2 ) omit what appears to be the governing 
physics of the spectral functions, namely rescattering and resonance effects. Consequently, the 
0(p 3 ) ChPT computation of the strange magnetic radius - though counterterm independent 
- contains only the non-analytic contributions at 0(g 2 ) and exceeds our upper bound for the 
magnetic radius by an order of magnitude. The higher-order rescattering corrections needed to 
render the ChPT prediction consistent with our bound are presumably contained in terms of 
0(p 4 ) or higher. Similarly, we suspect that the pattern of cancellations obtained in the 0(g 2 ) 
NRQM calculation of Ref. |9j will be significantly modified when rescattering and resonance 
effects are included. 

A computation to all orders in g of the remaining contributions would appear to be a daunting 
task. A few observations may simplify the problem, however. First, unitarity arguments suggest 
that the important structure in the spectral function lies below the two-nucleon threshold. 
Contributions from states such as AA, AA71", . . . whose thresholds t\ > Am 2 N are limited by 
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unitarity bounds on the strong partial wavesQ and are unlikely to be significantly enhanced by 
resonance effects in the intermediate state form factors. 

Second, we expect that the only important pionic contributions arise via resonances, such 
as the u;(780) or 0(1020). In ChPT, for example, the matrix element(37r|s7 At s|0) receives no 
non-analytic contributions at 0(p 7 ). Consequently, the 3ir contribution to the spectral function 
is small in the absence of resonant short distance effects. As discussed in Ref . , such effects 



(e.g., 37r <-> oj and 3tt <-> pir 0) may enhance the 3ir contribution up to the scale of the KK 
contribution. In fact, we speculate that the strength obtained in the VMD analyses of the 
isoscalar EM form factors arises primarily in the 3tc channel. There exists little evidence for 
significant coupling of higher-mass multi-pion states to [I G (J PC ) = 0~(1 )] resonances. We 
thus expect their contributions to be no larger than the non-resonant part of the KK term. 

Third, states involving pions and strange mesons may generate important contributions via 
the a; (1420), u;(1600) and 0(1680) resonances. A preliminary exploration of this possibility is 
given in Ref. |38| . In the VMD fits of Refs. |15|, [16], [I7| , inclusion of a vector meson pole in this 
mass region is needed to obtain an acceptable \ 2 . Since the flavor content of the vector mesons 
in this region is not known, the higher mass contributions to the strangeness form factors have 
been inferred from a priori assumptions about their large-t behavior. A reasonable range for 



the strange moments avoiding these assumptions has recently been given in Ref. |37 



A calculation of unitarity bounds for states having t\ > 4m N is tractable. Data for NN — > 3n 
and tcN — > ttttN may permit a model-independent determination of the 3tc contribution to all 
orders in g. Whether or not a realistic treatment of the other multi-meson states can be carried 
out remains to be seen. 
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